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[i]  We  investigate  the  solar  cycle  modulation  of  the  quasi-biennial  oscillation  (QBO) 
in  stratospheric  zonal  winds  and  its  impact  on  stratospheric  ozone  with  an  updated  version 
of  the  zonally  averaged  CHEM2D  middle  atmosphere  model.  We  find  that  the  duration 
of  the  westerly  QBO  phase  at  solar  maximum  is  3  months  shorter  than  at  solar 
minimum,  a  more  robust  result  than  in  an  earlier  CHEM2D  study  due  to  reduced  Rayleigh 
friction  drag  in  the  present  version  of  the  model.  The  modeled  solar  cycle  ozone 
response,  determined  via  multiple  linear  regression,  is  compared  with  observational 
estimates  from  the  combined  Solar  Backscattered  Ultraviolet  (SBUV/2)  data  set  for  the 
period  1979-2003.  We  find  that  a  model  simulation  including  imposed  solar  UV 
variations,  the  zonal  wind  QBO,  and  an  imposed  11 -year  variation  in  planetary 
wave  1  amplitude  produces  a  lower  stratospheric  ozone  response  of  ~2.5%  between 
0  and  20°S  and  an  upper  stratospheric  ozone  response  of  ~1%  between  45  and  55  km, 
in  good  agreement  with  the  SBUV-derived  ozone  response.  This  simulation  also 
produces  an  (enhancement/reduction)  in  the  (lower/upper)  stratospheric  temperature 
response  at  low  latitudes  compared  to  the  effects  of  solar  UV  variations  alone, 
which  are  consistent  with  model  vertical  velocity  anomalies  produced  by  the 
solar-modulated  QBO  and  imposed  changes  in  planetary  wave  forcing. 
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in  a  zonally  averaged  photochemical  transport  model  of  the  middle  atmosphere,  J.  Geophys.  Res.,  112,  D16109,  doi:10.1029/ 
2006JD008369. 


1.  Introduction 

[2]  A  long-standing  issue  in  understanding  Sun-climate 
connections  is  that  both  the  magnitude  and  the  altitude 
dependence  of  the  apparent  solar  cycle  variations  in  strato¬ 
spheric  ozone  derived  from  satellite  data  sets  differ  markedly 
from  estimates  based  on  current  models  of  stratospheric 
photochemistry  [. McCormack  and  Hood,  1996;  Shindell 
et  al.,  1999;  Hood,  2004],  Standard  photochemical  theory 
indicates  that  higher  levels  of  ultraviolet  (UV)  irradiance 
during  the  maximum  in  the  11 -year  solar  cycle  increase 
the  photolysis  rate  of  molecular  oxygen  to  produce  more 
stratospheric  ozone,  and  thus  more  radiative  heating,  in  the 
upper  stratosphere  (35-50  km  altitude)  through  ozone  UV 
absorption.  Photochemical  model  estimates  of  the  low- 
latitude  ozone  response  to  imposed  changes  in  solar  UV 
all  generally  peak  at  2%  near  35-40  km,  and  show 
variations  of  1%  or  less  both  above  and  below  this  region 
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\Hood,  2004;  Soukharev  and  Hood,  2006].  In  contrast, 
the  11 -year  ozone  response  derived  from  three  different 
satellite-based  data  sets  finds  a  different  pattern  in  the 
altitude  dependence  of  the  ozone  response  at  low  latitudes 
[ Soukharev  and  Hood,  2006].  This  pattern  consists  of  a 
positive  response  on  the  order  of  2-4%  (i.e.,  higher  ozone 
at  solar  maximum)  in  the  upper  stratosphere  between 
40  and  45  km,  a  zero  or  weakly  negative  response  in 
the  middle  stratosphere  between  30  and  40  km,  and  a 
positive  response  of  2-4%  in  the  lower  stratosphere 
between  20  and  30  km. 

[3]  Numerous  modeling  studies  have  shown  that  solar- 
induced  changes  in  stratospheric  ozone  heating  can  produce 
changes  in  temperatures  and  winds  that  alter  the  upward 
propagation  of  planetary-scale  waves  forced  in  the  extra- 
tropical  troposphere,  ultimately  producing  a  dynamical 
feedback  that  amplifies  the  original  solar-induced  effect 
and  translates  it  to  lower  levels  of  the  atmosphere  [e.g., 
Haigh,  1996;  Kodera  et  al.,  1997;  Rind,  2002;  Kodera  and 
Kuroda,  2002;  Matthes  et  al.,  2004,  2006],  Consistent  with 
this  hypothesis,  global  data  sets  of  stratospheric  ozone, 
temperature,  and  zonal  wind  exhibit  quasi-decadal  variability 
that  is  approximately  in  phase  with  the  11 -year  solar  cycle 
[e.g.,  Labitzke  and  van  Loon,  1988;  Hood  et  al.,  1993; 
Kuroda  and  Kodera,  2002;  World  Meteorological  Organi- 
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zation  ( WMO ),  2003;  Hood,  2004;  Crooks  and  Gray,  2005], 
More  recently,  Austin  et  al.  [2006]  report  that  variations 
in  sea  surface  temperatures  may  also  contribute  to  the 
observed  altitude  dependence  of  the  ozone  response.  The 
discrepancy  between  the  predicted  and  observed  ozone 
response  to  solar  UV  variability  is  a  major  source  of 
uncertainty  in  our  understanding  of  possible  mechanisms 
linking  solar  variability  and  climate.  The  present  study  uses 
a  zonally  averaged  (2D)  photochemical  transport  model 
of  the  middle  atmosphere  to  identify  possible  dynamical 
origins  of  this  discrepancy. 

[4]  Currently,  it  is  unclear  how  much  of  the  observed 
altitude  dependence  in  the  ozone  solar  response  is  a  true 
geophysical  response  to  solar  UV  changes  [e.g.,  Hood, 
2004],  and  how  much  may  be  a  statistical  artifact  arising 
from  the  relatively  short  length  of  existing  global  satellite 
data  sets  and  the  occurrence  of  large  volcanic  eruptions  near 
periods  of  solar  maximum  [e.g.,  Lee  and  Smith,  2003], 
Since  stratospheric  ozone  photochemistry  is  relatively  well 
understood,  dynamical  processes  not  included  in  most 
photochemical  models  are  often  cited  as  a  possible  geo¬ 
physical  origin  for  the  observed  vertical  structure  in  the 
ozone  response.  Earlier  studies  have  shown  that  the  positive 
lower  stratospheric  ozone  response  is  separate  and  distinct 
from  volcanic  effects  [ McCormack  et  al.,  1997]  and  is 
consistent  with  a  quasi-decadal  modulation  in  the  dynamics 
of  the  tropical  lower  stratosphere  [Hood,  1997;  Hood  and 
Soukharev,  2003].  In  the  present  study,  we  explore  the 
effects  of  solar-induced  dynamical  variations  on  lower 
stratospheric  ozone  with  the  zonally  averaged  CHEM2D 
photochemical  transport  model. 

[5]  Support  for  a  solar  influence  on  lower  stratosphere 
dynamics  is  found  in  the  observed  11 -year  variation  of  the 
quasi-biennial  oscillation  (QBO)  in  stratospheric  equatorial 
zonal  winds.  A  quasi-decadal  modulation  of  the  westerly 
QBO  phase  was  first  reported  by  Quiroz  [1981]  and  later  by 
Salby  and  Callaghan  [2000]  using  an  extended  data  record. 
An  analysis  of  equatorial  zonal  winds  by  Hamilton  [2002] 
confirms  this  result,  but  concludes  that  available  data  sets 
are  still  of  insufficient  length  to  unambiguously  attribute 
this  variability  to  changes  in  solar  activity.  Pascoe  et  al. 
[2005]  report  evidence  of  a  solar  cycle  modulation  of  the 
mean  descent  rate  of  the  easterly  shear  zone  between  20  and 
44  hPa  in  the  ERA-40  data  set,  such  that  the  easterly  shear 
descends  about  2  month  faster  during  solar  maximum. 
This  is  qualitatively  consistent  with  the  shorter  westerly 
QBO  phase  during  solar  maximum  reported  by  Salby  and 
Callaghan  [2000]. 

[6]  There  also  exists  a  complex  relationship  between  the 
QBO  and  planetary  wave  activity  in  the  wintertime  extra- 
tropical  stratosphere  that  may  also  be  modified  by  the  solar 
cycle.  For  example,  breaking  planetary  waves  in  the 
wintertime  stratosphere  drive  the  mean  meridional  circula¬ 
tion  in  the  middle  atmosphere  [ Holton  et  al,  1995],  and  the 
propagation  of  these  waves  from  the  troposphere  is  sensitive 
to  the  phase  of  the  zonal  wind  QBO  [Baldwin  et  al.,  2001]. 
In  addition,  the  descent  of  the  QBO  easterlies  is  known  to 
occur  preferentially  during  the  months  of  May -July,  when 
the  vertical  component  of  the  mean  meridional  circulation  in 
the  tropical  lower  stratosphere  is  weakest.  A  modeling  study 
by  Gray  et  al.  [2004]  suggests  that  equatorial  upper  strato¬ 
spheric  zonal  wind  anomalies  related  to  both  the  solar  cycle 


and  the  QBO  can  alter  planetary  wave  propagation  and  thus 
affect  the  frequency  and  timing  of  stratospheric  sudden 
warmings  in  the  Northern  Hemisphere  (NH)  winter.  The 
present  study  offers  evidence  that  the  combined  effects 
of  solar  UV  variations,  the  QBO,  and  planetary  wave 
forcing  could  provide  a  possible  explanation  for  the 
observed  11 -year  ozone  response  in  the  tropical  lower 
stratosphere,  where  the  ozone  distribution  is  largely 
controlled  by  transport. 

[7]  2D  modeling  studies  are  useful  for  evaluating  the 
possible  solar  cycle  modulation  of  the  QBO,  since  the 
relative  simplicity  of  the  2D  model  makes  it  straightforward 
to  isolate  internal  and  external  sources  of  interannual 
variability  over  multidecadal  periods.  Similar  calculations 
with  a  3D  general  circulation  model  (GCM)  are  much  more 
computationally  expensive,  and  the  results  are  more  diffi¬ 
cult  to  interpret.  For  example,  while  a  recent  GCM  study  by 
Palmer  and  Gray  [2005]  using  an  internally  generated  zonal 
wind  QBO  found  the  QBO  period  to  be  2  months  shorter  at 
solar  maximum  than  at  solar  minimum,  an  unrealistically 
large  change  in  solar  irradiance  was  needed  to  produce  a 
statistically  significant  response  in  the  relatively  short 
model  integration  (25  years). 

[8]  The  2D  modeling  study  of  McCormack  [2003],  using 
an  interactive  representation  of  the  zonal  wind  QBO,  first 
showed  that  realistic  changes  in  solar  UV  irradiance  are 
capable  of  modulating  the  QBO  period  in  a  manner  that  was 
qualitatively  consistent  with  the  observed  effect,  i.e.,  a 
shorter  duration  of  the  QBO  westerly  phase  during  solar 
maximum.  However,  the  difference  in  the  duration  of  the 
modeled  QBO  westerly  phase  between  solar  minimum  and 
solar  maximum  was  only  1  month,  whereas  observations 
indicate  a  difference  of  3-6  months  [Salbv  and  Callaghan, 
2000], 

[9]  The  present  study  uses  an  updated  version  of  the 
CHEM2D  model  used  by  McCormack  [2003]  to  simulate 
the  interannual  variability  in  stratospheric  ozone  in  response 
to  the  11 -year  solar  cycle,  the  QBO,  and  quasi-decadal 
changes  in  the  dynamical  forcing  of  the  stratospheric 
meridional  circulation  by  planetary  waves.  This  version  of 
CHEM2D  features  a  higher  model  top,  increased  vertical 
resolution,  and  an  improved  treatment  of  gravity  wave  drag 
in  the  middle  atmosphere  in  which  an  imposed  Rayleigh- 
type  drag  on  the  zonal  winds  now  plays  a  lesser  role 
[Siskind  et  al.,  2003].  The  goal  of  this  study  is  to  determine 
if  the  observed  vertical  structure  in  the  11 -year  ozone 
response  can  be  explained  by  the  interaction  of  the  solar 
cycle,  QBO,  and  planetary  wave  effects.  We  will  focus  in 
particular  on  the  role  of  the  QBO  in  producing  the  secondary 
maximum  in  the  observed  11 -year  ozone  response  in  the 
tropical  lower  stratosphere  (i.e.,  15-25  km  altitude).  Our 
results  indicate  that  the  improved  model  formulation  produ¬ 
ces  a  more  robust  solar  cycle  modulation  of  the  zonal  wind 
QBO  compared  to  earlier  results  [McCormack,  2003],  We 
also  find  that  the  combined  effects  of  solar  UV  variations,  the 
zonal  QBO,  and  planetary  wave  forcing  can  enhance  the 
ozone  response  to  changes  in  solar  UV  irradiance  in  both 
the  lower  and  upper  stratosphere  in  a  manner  that  is  qualita¬ 
tively  consistent  with  observations. 

[10]  A  description  of  the  CHEM2D  model  is  given  in 
section  2.  Section  3  describes  the  analysis  of  model  output 
using  multiple  linear  regression.  Section  4  presents  results 
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Figure  1.  (a)  Altitude/time  section  of  CHEM2D  equatorial  zonal  wind  from  the  control  experiment 
(EXPO),  solid  (dashed)  contours  denoting  westerly  (easterly)  winds  are  drawn  for  0,  ±10,  ±20,  ±40,  and 
±60  m  s  ',  (b)  Vertical  profiles  of  annual  mean  CHEM2D  equatorial  zonal  wind  (solid  curve)  and 
climatological  values  obtained  from  the  UARS  Reference  Atmosphere  Project  (dashed  curve). 


from  a  series  of  CHEM2D  simulations  including  the  effects 
of  the  11 -year  solar  cycle  and  the  QBO.  Section  5  discusses 
how  these  effects  produce  better  agreement  between  the 
observed  and  modeled  solar  cycle  variations  in  ozone. 
Section  6  contains  a  summary  and  outlines  possible  avenues 
of  future  research. 

2.  CHEM2D  Model 

2.1.  Model  Description 

[n]  CHEM2D  is  a  zonally  averaged  (2D)  model  that 
features  a  fully  self-consistent  treatment  of  radiative,  pho¬ 
tochemical,  and  dynamical  processes  of  the  middle  atmo¬ 
sphere.  [McCormack  and  Siskind,  2002;  Siskind  et  al, 
2003;  McCormack  et  al,  2006].  The  model  photochemistry 
accounts  for  54  different  species,  and  now  includes  catalytic 
ozone  loss  due  to  bromine  compounds.  Model  reaction  rates 
are  based  on  Sander  et  al.  [2003],  The  radiative  and 
photochemical  calculations  are  performed  once  per  day, 
and  the  model  dynamics  are  updated  every  2  hours. 
The  model  extends  from  pole  to  pole  with  grid  points 
spaced  every  4.8°  in  latitude.  The  model  vertical  domain 
consists  of  88  levels  from  the  surface  to  pressure  level 
p  —  6  x  10' 5  hPa  (~116  km)  that  are  spaced  every 
1.3  km.  This  version  of  the  model  offers  twice  the  vertical 
resolution  as  earlier  versions  [McCormack  and  Siskind, 
2002;  McCormack,  2003]. 

[12]  The  CHEM2D  dynamical  framework  is  based  on  the 
Transformed  Eulerian  Mean  formulation,  in  which  the 
residual  meridional  circulation  is  driven  by  zonally  aver¬ 
aged  sources  of  momentum  and  thermodynamic  forcing 
[McCormack  and  Siskind,  2002],  The  momentum  sources 
include  a  parameterization  for  subgrid-scale  gravity  wave 
drag  from  both  stationary  and  nonzero  phase  speed  gravity 
waves,  the  Eliassen-Palm  (EP)  flux  divergence  associated 
with  dissipating  or  breaking  planetary  waves,  and  parame¬ 
terized  vertical  mixing  of  momentum  by  gravity  wave 
breaking  and  molecular  diffusion  in  the  upper  mesosphere. 
An  altitude-dependent  Rayleigh  friction  term  is  used  to 
constrain  the  zonal  wind  field  near  the  model  top. 


[13]  To  simulate  the  zonal  wind  QBO,  the  CHEM2D 
model  includes  a  parameterization  of  the  upward  momen¬ 
tum  flux  divergence  associated  with  equatorial  Kelvin  and 
Rossby-gravity  wave  modes  [Holton  and  Lindzen,  1972; 
Dunkerton,  1979;  Gray  and  Pyle,  1989].  For  a  detailed 
description  of  the  CFIEM2D  QBO  parameterization,  see 
McCormack  and  Siskind  [2002],  This  parameterization 
allows  the  period  of  the  QBO  to  vary  in  response  to  changes 
in  the  strength  of  the  model  residual  circulation.  It  has  been 
used  previously  in  the  CEIEM2D  model  to  simulate  QBO- 
related  features  in  stratospheric  constituent  transport 
[McCormack  and  Siskind,  2002]  and  interactions  between 
solar  UV  irradiance  changes  and  the  period  of  the  zonal 
wind  QBO  [McCormack,  2003],  Figure  la  plots  an  altitude¬ 
time  section  of  CHEM2D  equatorial  zonal  winds  illustrating 
the  characteristics  of  the  modeled  QBO.  The  period  of  the 
QBO  ranges  from  24  to  33  months  with  an  average  period 
of  28  months,  similar  to  observations. 

[14]  One  advantage  of  the  CFIEM2D  model’s  interactive 
QBO  parameterization  is  that,  unlike  models  that  employ  a 
sinusoidal  momentum  source  with  a  fixed  period  [e.g.,  Lee 
and  Smith,  2003],  it  can  capture  possible  solar-related 
photochemical-dynamical  feedbacks  that  may  influence 
the  QBO.  Models  that  produce  a  zonal  wind  QBO  by 
relaxation  to  observed  equatorial  winds  [e.g.,  Matthes  et 
al,  2004]  may  capture  the  solar-QBO  interaction,  but  they 
cannot  isolate  and  quantify  the  origin  of  the  interaction. 

[15]  Another  advantage  of  the  CHEM2D  model  is  that  it 
is  relatively  inexpensive  for  performing  long-term  climate 
simulations.  Recent  GCM  studies  have  been  able  to  produce 
realistic  QBO  behavior  in  model  zonal  winds  either  by 
explicitly  resolving  equatorial  wave  modes  with  high  ver¬ 
tical  resolution  [Giorgetta  et  al,  2006]  or  by  adjusting  the 
subgrid-scale  gravity  wave  drag  parameterization  [Palmer 
and  Gray,  2005],  To  fully  investigate  the  impact  of  solar- 
QBO  interactions  on  stratospheric  ozone,  these  type  of 
simulations  should  extend  50  years  or  longer  and  include 
fully  coupled  photochemistry,  realistic  solar  UV  forcing, 
and  a  realistic  zonal  wind  QBO. 
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Table  1.  Description  of  Experiments 


Experiment 

Description 

EXPO 

fixed  solar  UV  with  zonal  wind  QBO 

EXP  I 

varying  solar  UV  without  zonal  wind  QBO 

EXP1RF 

same  as  EXP1  using  McCormack  [2003] 
Rayleigh  friction  profile 

EXP2 

varying  solar  UV  with  zonal  wind  QBO 

EXP3 

varying  solar  UV  with  zonal  wind  QBO  and 
planetary  wave  forcing 

EXP4 

varying  solar  UV  with  planetary  wave 
forcing  without  QBO 

[i6]  The  CHEM2D  gravity  wave  drag  parameterization 
includes  both  stationary  waves  and  nonstationary  waves 
with  discrete  phase  speeds  ranging  from  —  50m/sto+50m/s 
having  a  specified  source  distribution  that  varies  with 
latitude  [ Siskind  et  al.,  2003].  Additional  nonstationary 
wave  modes  have  been  introduced  in  the  tropics  to  produce 
a  realistic  semiannual  oscillation  (SAO)  in  equatorial  zonal 
wind  throughout  the  middle  atmosphere  [ McCormack , 
2003],  Figure  lb  compares  the  annual  average  CHEM2D 
zonal  wind  profile  over  the  equator  with  a  climatological 
wind  profile  based  on  observations  from  the  Upper  Atmo¬ 
sphere  Research  Satellite  (UARS).  In  general,  the  CHEM2D 
zonal  winds  show  good  agreement  with  the  UARS  clima¬ 
tology,  although  there  is  an  easterly  bias  in  the  CHEM2D 
equatorial  zonal  winds  between  40  and  50  km  altitude. 

[n]  The  increased  model  vertical  resolution  in  the  present 
version  of  CHEM2D  was  motivated  in  part  by  the  study  of 
McCormack  [2003],  which  found  a  solar  cycle  modulation 
of  the  zonal  wind  QBO  period  that  was  qualitatively 
consistent  with,  but  smaller  than,  the  quasi-decadal  modu¬ 
lation  reported  by  Salby  and  Callaghan  [2000].  To  improve 
the  agreement  between  the  observed  and  modeled  solar- 
QBO  effect,  we  have  increased  the  model  vertical  resolution 
to  obtain  a  better  representation  of  parameterized  subgrid- 
scale  wave  drag  on  the  zonal  winds.  Results  presented  in 
section  4.2  show  that  by  decreasing  the  effect  of  Rayleigh 
friction  on  zonal  winds  near  the  stratopause,  the  present 
high-vertical-resolution  CHEM2D  model  produces  a  more 
robust  solar  cycle  modulation  of  the  QBO  period  compared 
to  the  earlier  result  of  McCormack  [2003]. 

[is]  CHEM2D  solar  UV  variations  over  the  1 1-year  cycle 
are  specified  as  a  function  of  wavelength  from  125  to 
300  nm  using  the  satellite-based  estimates  of  Lean  et  al 
[1997].  These  variations  are  imposed  as  a  sinusoidal  1 1 -year 
variation.  Unlike  earlier  studies  that  have  performed  steady 
state  calculations  under  both  solar  maximum  and  solar 
minimum  conditions,  this  time-varying  approach  allows  us 
to  evaluate  the  solar  response  through  multiple  realizations 
of  the  QBO  and  annual  cycle.  This  is  important  since 
the  phase  relationship  between  the  QBO  and  annual 
cycle  changes  over  time,  which  can  introduce  an  internal 
beat  period  in  model  transport  fields  on  decadal  timescales 
[e.g.,  McCormack  and  Siskind,  2002,  Figure  11].  Further¬ 
more,  a  recent  3D  modeling  study  by  Austin  et  al.  [2006] 
demonstrates  that  using  time-varying  solar  forcing  produces 
a  stratospheric  ozone  response  that  is  in  better  agreement 
with  observations  than  using  fixed  solar  forcing  for 


maximum  and  minimum  conditions  separately,  even  when 
the  effects  of  the  QBO  are  not  considered. 

2.2.  Experimental  Setup 

[19]  Table  1  lists  a  series  of  CFIEM2D  model  experiments 
that  are  designed  to  further  examine  the  solar  cycle  modu¬ 
lation  of  the  zonal  wind  QBO  reported  by  McCormack 
[2003]  and  to  assess  the  combined  effects  of  solar  UV 
variations,  the  QBO,  and  planetary  wave  activity  on  the 
11 -year  response  in  stratospheric  ozone  and  temperature. 
All  model  simulations  are  50  years  in  length  unless 
otherwise  noted. 

[20]  We  begin  with  a  control  experiment  (EXPO  in 
Table  1)  that  includes  the  interactive  QBO  parameteriza¬ 
tion  with  fixed  solar  UV  irradiance.  In  the  first  experiment 
with  varying  solar  UV  (EXP1),  we  impose  an  11 -year 
sinusoidal  variation  in  solar  UV  irradiance  based  on 
estimates  from  Lean  et  al.  [1997].  The  zonal  wind  QBO 
parameterization  is  not  included  in  EXP1  in  order  to 
produce  a  “pure”  solar  signal  in  the  model  ozone  fields. 
An  additional  simulation  without  the  QBO  (EXPIRE)  uses 
the  Rayleigh  friction  drag  profile  of  McCormack  [2003]  to 
demonstrate  how  it  affects  the  model’s  dynamical  response 
to  solar  UV  variations  in  the  upper  stratosphere  (see 
section  4.1). 

[21]  Experiment  2  (EXP2)  is  an  extended  150-year  sim¬ 
ulation  that  includes  both  the  11 -year  solar  UV  variations 
and  the  zonal  wind  QBO  parameterization.  With  the  longer 
simulation,  we  are  able  to  obtain  a  statistically  significant 
solar  cycle  modulation  of  the  QBO  (see  section  4.2). 
Experiment  3  (EXP3)  introduces  an  11 -year  modulation  in 
the  extratropical  planetary  wave  forcing,  which  is  consistent 
with  the  observational  estimates  of  Hood  and  Soukharev 
[2003],  in  addition  to  the  solar  UV  and  QBO  effects.  Finally, 
the  fourth  experiment  (EXP4)  includes  the  solar  UV  and 
planetary  wave  forcing  effects  but  no  zonal  wind  QBO 
parameterization.  A  comparison  of  results  from  EXP3  and 
EXP4  shows  that  the  imposed  11 -year  variations  in  solar 
UV  and  planetary  wave  amplitudes  produce  a  larger  quasi- 
decadal  signal  in  lower  stratospheric  ozone  when  the  zonal 
wind  QBO  is  present. 

3.  Linear  Regression  Model 

[22]  This  section  describes  the  multiple  regression  statis¬ 
tical  model  used  to  analyze  the  results  from  the  experiments 
listed  in  Table  1.  Multiple  linear  regression  models  have 
long  been  a  standard  tool  for  estimating  trends  in  strato¬ 
spheric  ozone  [e.g.,  WMO,  1991;  Harris  et  al.,  1998,  and 
references  therein].  The  standard  approach  is  to  assume 
that  the  temporal  behavior  of  a  zonally  averaged  time 
series  (i.e.,  ozone  or  temperature)  can  be  approximated  by 
a  least-squares  fit  to  a  linear  model  with  trend,  QBO,  and 
11 -year  solar  cycle  components.  Initially,  the  motivation 
for  this  approach  was  to  ensure  that  natural  sources  of 
interannual  variability  did  not  project  onto  estimates  of  the 
anthropogenic  trends  in  ozone.  With  the  advent  of  longer 
satellite -based  data  sets,  it  became  feasible  to  apply 
this  model  to  quantify  the  effects  of  the  QBO  and  the  solar 
cycle  on  stratospheric  ozone  and  temperature  in  order  to 
facilitate  model-data  comparisons  of  interannual  variability 
[Hood  and  McCormack,  1992;  Randel  and  Cobb,  1994; 
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Year 

Figure  2.  Standardized  time  series  representing  solar  UV 
variations  and  the  first  two  principal  components  of  the 
zonal  wind  QBO  used  in  the  multiple  linear  regression 
model. 

McCormack  and  Hood,  1996;  Crooks  and  Gray,  2005; 
Soukharev  and  Hood,  2006], 

[23]  Here  we  adopt  an  approach  similar  to  that  of  Lee  and 
Smith  [2003]  where  we  apply  a  multiple  linear  regression 
model  to  output  from  a  2D  model  in  order  to  isolate  the 
11 -year  variability  in  ozone  and  temperature.  The  regression 
model  quantifies  the  components  of  interamiual  variability 
in  the  dependent  time  series  (i.e.,  the  model  output)  asso¬ 
ciated  with  a  set  of  independent  explanatory  time  series 
(i.e.,  the  11 -year  solar  UV  forcing  and  the  QBO).  Using 
this  approach,  we  are  able  to  make  direct  comparisons 
between  the  modeled  solar  cycle  variability  in  ozone 
and  the  SBUV(/2)-derived  estimates  of  the  11 -year  ozone 
response  of  Soukharev  and  Hood  [2006]. 


[24]  One  difference  between  our  regression  model  and 
that  of  Soukharev  and  Hood  [2006]  is  the  absence  of  a  trend 
term  and  the  accompanying  autocorrelation  term.  Aside 
from  the  solar  UV  and  QBO  effects  described  above,  no 
interannual  variability  or  trends  are  directly  imposed  on  the 
model’s  photochemical  constituents,  therefore  there  is  no 
need  to  model  this  trend  statistically.  Regression  analyses  of 
CHEM2D  model  ozone  output  conducted  both  with  and 
without  the  trend  and  autocorrelation  terms  show  negligible 
differences  in  the  resulting  estimates  of  the  solar  cycle 
response.  Therefore  we  elect  not  to  include  this  term  in 
our  analysis. 

[25]  Previous  observational  studies  have  described  the 
QBO  in  their  regression  models  using  either  a  lagged  time 
series  of  the  zonal  wind  QBO  at  one  particular  level 
[McCormack  and  Hood,  1996;  Soukharev  and  Hood, 
2006],  a  combination  of  zonal  wind  time  series  at  two 
different  levels  where  the  time  variability  is  in  approximate 
quadrature  [Brunner  et  al.,  2006],  or  proxy  time  series  based 
on  the  set  of  linear  orthogonal  basis  functions  determined 
from  a  principal  component  or  empirical  orthogonal  func¬ 
tion  analysis  of  equatorial  stratospheric  zonal  winds  [e.g., 
Randel  and  Wu,  1996;  Crooks  and  Gray,  2005],  We  have 
tested  each  form  of  QBO  representation  in  our  regression 
model  and  found  no  significant  differences  in  the  derived 
solar  regression  coefficients.  We  note  that  the  CHEM2D 
model’s  representation  of  the  zonal  wind  QBO  is  only  an 
approximation  of  the  true  zonal  wind  QBO,  and  as  such 
may  not  capture  all  the  temporal  variability  in  the  observed 
QBO.  Recent  2D  photochemical  modeling  studies  have 
shown  that  a  lagged  QBO  proxy  time  series  may  introduce 
artifacts  in  solar  regression  coefficients  if  the  QBO  period 
is  highly  irregular  over  the  course  of  several  decades 
(A.  Smith,  personal  communication,  2006).  To  avoid  this 
possibility,  our  regression  model  employs  a  QBO  proxy 
time  series  based  on  the  results  of  a  principal  component 
analysis  of  CHEM2D  equatorial  zonal  winds. 


Latitude  Latitude 


Figure  3.  Spatial  patterns  of  the  ozone  regression  coefficients  (a)  Jqboi  and  (b)  Jqbo2,  associated  with 
the  first  and  second  principal  components  of  the  CHEM2D  equatorial  stratospheric  zonal  wind, 
respectively,  from  the  control  experiment  (EXPO).  Coefficients  are  expressed  as  the  percent  change  for  a 
one  standard  deviation  change  in  the  principal  component  time  series.  Shaded  regions  are  not  statistically 
significant  at  the  2a  level. 
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CHEM2D  03  SOLAR  COEF  (%,  Max  -  Min) 


Figure  4.  Annual  average  ozone  solar  regression  coeffi¬ 
cients,  expressed  as  the  percent  change  (solar  maximum 
minus  solar  minimum)  from  experiment  1,  which  includes 
an  imposed  1 1 -year  cycle  in  solar  UV.  Solid  contours  denote 
positive  values,  dashed  contours  denote  negative  values, 
and  shaded  areas  are  not  significant  at  the  2a  level. 
Contours  are  drawn  every  ±0.5%  up  to  ±4%,  and  then  every 
±2%  thereafter. 

[26]  The  multiple  linear  regression  model  used  in  the 
present  study  is  of  the  form 

y(t)  =  A/(0  +  lsxs{t)  +  ~1qbo\xQbo\  +  1 IqboiXqboi •  (1) 

Here  y  is  the  model  time  series  (i.e.,  zonal  mean  ozone  or 
temperature),  /(,  is  the  corresponding  long-term  mean  value 
for  each  month  i  =  1,. . .,  12,  and  Xs  is  a  11-year  sinusoidal 
variation  in  solar  UV  (see  Figure  2,  top  curve).  XqBO\  and 
Xqboi  are  time  series  representing  the  two  leading  principal 
components  derived  from  a  linear  system  of  model 
equatorial  zonal  winds  between  10  and  70  hPa  from  each 
experiment.  The  ozone  (temperature)  regression  coefficients 
7s,  7 qbou  and  7gso2  are  calculated  in  terms  of  volume 
mixing  ratio  (Kelvin)  on  a  constant  pressure  (p )  grid 
and  converted  to  altitude  (z)  using  the  relation  z  = 
— F/log(/V1013),  where  the  pressure  scale  height  H  is 
taken  to  be  a  constant  7.021  km.  The  seasonal  dependence 
of  each  regression  coefficient  is  represented  by  a  series  of 
harmonic  functions  [Randel  and  Cobb,  1994].  Uncertainties 
in  the  regression  coefficients  are  determined  from  the 
error  covariances  of  the  least-squares  fit.  Where  appropriate, 
±2cr  values  of  a  regression  coefficient’s  uncertainties  are 
quoted  in  parentheses. 

[27]  To  illustrate  the  QBO  proxy  time  series,  Figure  2 
plots  the  normalized  time  series  of  XqBO  1  and  XqB02  from 
the  control  run  (EXPO).  Together  they  explain  over  90%  of 
the  temporal  variance  in  the  model  equatorial  zonal  winds 
between  10  and  70  hPa.  Figure  3  shows  the  latitude  and 
altitude  dependences  of  the  ozone  regression  coefficients 
1qbo\  and  7 qBoi  calculated  from  regression  analysis  of  the 
CHEM2D  ozone  fields  from  the  EXPO  run;  shaded  regions 
are  not  statistically  significant  at  the  2 a  level.  For  each 
experiment,  the  QBO  proxy  time  series  are  determined  by 


principal  component  analysis  of  output  zonal  wind  fields 
prior  to  the  regression  analysis. 

[28]  The  regression  model  (1)  is  applied  to  the  CHEM2D 
model  ozone  and  temperature  fields  from  each  of  the 
experiments  listed  in  Table  1,  which  are  diumally  averaged 
and  output  on  the  first  day  of  each  calendar  month.  All 
regression  analyses  are  based  on  50  years  of  model  output 
to  facilitate  comparison  of  regression  coefficients’  statistical 
significance  among  the  different  experiments. 

4.  Results 

[29]  Here  we  present  results  from  the  CHEM2D  model 
experiments  that  are  described  in  section  2.2  and  summa¬ 
rized  in  Table  1. 

4.1.  Experiment  1:  Varying  Solar  UV  Only 

[30]  Figure  4  plots  the  annual  average  solar  regression 
coefficients  derived  from  the  EXP1  ozone  fields  between 
1 8  and  65  km,  which  are  expressed  in  terms  of  the  percentage 
change  in  ozone  from  solar  minimum  to  maximum.  Consis¬ 
tent  with  earlier  modeling  studies  [Hood,  2004],  we  find 
the  largest  positive  ozone  response  (~3%)  between  35  and 
40  km,  in  contrast  to  the  observational  results  that  find  the 
largest  response  in  the  upper  stratosphere.  Shaded  regions  in 
Figure  4  denote  where  the  solar  regression  coefficients  from 
this  50-year  run  are  not  different  from  zero  at  the  2a  level. 
Typical  2a  uncertainty  values  from  the  analysis  of  EXP  1  are 
0.1%  or  less.  As  Figure  4  indicates,  the  EXP1  regression 
coefficients  are  statistically  significant  throughout  most  of 
the  stratosphere  and  lower  mesosphere. 

[31]  Next,  we  compare  the  modeled  stratospheric  ozone 
response  with  observational  estimates  derived  from 
the  combined  24-year  (1979-2003)  Solar  Backscatter 
Ultraviolet  (SBUV)  and  SBUV2  ozone  profile  data  set 
by  Soukharev  and  Hood  [2006]  using  a  similar  multiple 


SBUV/2  A  03  (%,  max-min)  :  1979-2003 


Figure  5.  Annual  average  ozone  solar  regression  coeffi¬ 
cients  derived  from  the  combined  SBUV(/2)  ozone  profile 
data  set  for  the  period  1979-2003  [Soukharev  and  Hood, 
2006],  Regression  coefficients  are  expressed  as  the 
percentage  difference  between  solar  maximum  and  solar 
minimum.  Shaded  regions  are  not  statistically  significant  at 
the  2a  level. 
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Figure  6.  (a)  Profiles  of  the  Rayleigh  friction  timescales  in  experiment  1  (black  curve)  and  experiment 
1RF  (gray  curve).  Experiment  1RF  is  identical  to  the  profile  used  by  McCormack  [2003],  (b)  Solar- 
induced  anomalies  in  the  equatorial  zonal  wind  at  50  km  from  experiment  1  (black  curve)  and  experiment 
1RF  (gray  curve).  Dashed  vertical  line  indicates  solar  maximum. 


linear  regression  model.  Figure  5  plots  the  annual  average 
SBUV(/2)  solar  regression  coefficients  as  a  function  of 
latitude  and  approximate  altitude.  Since  SBUV  profile 
measurements  are  limited  to  sunlit  regions,  solar  regression 
coefficients  are  only  computed  for  the  latitude  region  from 
60°S  to  60°N.  Typical  values  of  the  2cr  uncertainties  in  the 
SBUV-derived  solar  regression  coefficients  are  ±1% 
throughout  most  of  the  stratosphere. 

[32]  As  mentioned  in  the  Introduction,  the  altitude 
and  latitude  dependences  of  the  modeled  solar  response 
in  stratospheric  ozone  from  EXP1  in  Figure  4  show  marked 
disagreement  with  the  SBUV(/2)-derived  response  in 
Figure  5.  Most  noticeably,  the  modeled  ozone  response  in 
the  tropics  peaks  near  35  km,  whereas  the  SBUV(/2)  solar 
coefficients  in  the  tropics  peak  between  45  and  50  km. 
Furthermore,  SBUV(/2)  data  show  no  statistically  signifi¬ 
cant  ozone  response  over  the  solar  cycle  in  the  tropics  from 
30  to  45  km.  Two  regions  with  somewhat  reasonable 
agreement  are  found  near  35  km  between  40-50°S  and 
50-60°N  latitude,  where  the  EXP1  results  indicate  a  2-3% 
ozone  increase  from  solar  minimum  to  solar  maximum.  In 
the  lower  stratosphere  (18-25  km),  the  EXP1  results  show 
little  to  no  significant  ozone  response  throughout  the 
tropics,  whereas  the  SBUV(/2)  results  indicate  a  secondary 
subtropical  maximum  in  the  ozone  response  with  peak 
values  of  2.5-3%  at  the  lowest  analyzed  level  at  50  hPa 
(~21  km). 

[33]  Currently,  it  is  not  clear  why  the  observed  upper 
stratospheric  ozone  response  exceeds  photochemical  model 
estimates.  Langematz  et  al.  [2005]  proposed  that  this  feature 
could  be  a  response  to  solar-modulated  production  and 
transport  of  mesospheric  odd  nitrogen  to  the  tropical  stra- 
topause  region.  However,  an  analysis  of  HALOE  NOx 
records  over  the  last  solar  cycle  by  Hood  and  Soukharev 
[2006]  shows  no  significant  11 -year  variation  in  the  upper 
stratosphere  to  support  this  assertion.  As  shown  in  section  4.3, 
model  simulations  that  include  imposed  variations  in  plane¬ 
tary  wave  forcing  (EXP3  and  EXP4)  produce  a  somewhat 
larger  upper  stratospheric  ozone  response  than  solar  forcing 


alone  (EXP1),  in  better  agreement  with  the  observed  upper 
stratospheric  ozone  response  (Figure  5). 

[34]  The  lack  of  a  significant  observed  ozone  response  in 
the  middle  stratosphere  (i.e.,  30-40  km  altitude)  could  not 
be  reproduced  in  EXP1.  Lee  and  Smith  [2003]  have 
suggested  that  the  zero  or  slightly  negative  midstratospheric 
ozone  response  seen  in  the  observations  is  due  to  the  effects 
of  major  volcanic  eruptions  that  occurred  near  the  maxi¬ 
mum  of  solar  cycles  21  and  22.  Their  assertion  was  based 
on  a  zonally  averaged  (2D)  photochemical  model  simula¬ 
tion  in  which  the  anomalous  heating  due  to  volcanic 
eruptions  produces  a  quasi-decadal  variation  of  ^50%  in 
the  strength  of  the  westerly  zonal  wind  QBO  phase, 
contrary  to  observations.  Since  the  CHEM2D  model  does 
not  include  heterogeneous  chemical  reactions  on  strato¬ 
spheric  aerosols,  we  cannot  test  this  claim  directly.  We  note 
that  an  intercomparison  of  the  observed  11 -year  ozone 
response  in  SBUV(/2),  HALOE,  and  SAGEI/I1  data  by 
Soukharev  and  Hood  [2006]  finds  a  minimum  in  the 
stratospheric  ozone  response  near  35  km  in  all  three  data 
sets  regardless  of  whether  or  not  periods  of  high  volcanic 
aerosol  loading  were  included  in  the  analyses. 

[35]  One  major  difference  between  the  high-vertical- 
resolution  model  used  in  the  present  study  and  the  lower- 
vertical-resolution  model  used  in  McCormack  [2003]  is  in 
the  magnitude  of  the  prescribed  Rayleigh  friction  term  used 
to  constrain  the  zonal  winds  in  the  mesosphere.  With  the 
higher  model  top  and  increased  vertical  resolution,  it  is  no 
longer  necessary  to  use  such  a  strong  damping  factor 
on  mesospheric  winds.  Figure  6a  compares  the  Rayleigh 
friction  profile  used  by  McCormack  [2003]  (gray  curve) 
with  the  profile  used  in  the  present  study  (black  curve).  The 
old  Rayleigh  friction  profile  makes  a  nonzero  contribution 
to  the  momentum  budget  near  the  stratopause,  where  the 
solar-induced  heating  anomalies  are  largest.  In  contrast,  the 
effects  of  the  new  Rayleigh  friction  profile  are  limited  to 
the  upper  portion  of  the  mesosphere. 

[36]  The  consequences  of  the  reduced  Rayleigh  friction 
profile  on  the  solar-induced  zonal  wind  anomalies  near  the 
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Figure  7.  (a)  Time  series  of  CHEM2D  equatorial  zonal 
winds  at  32  hPa  from  experiment  1  with  fixed  solar  UV 
(black  curve)  and  experiment  2  with  varying  solar  UV  (gray 
curve).  Dashed  vertical  lines  indicate  solar  maximum. 

stratopause  can  be  seen  in  Figure  6b,  which  compares  the 
solar  cycle  response  in  equatorial  zonal  winds  at  50  km 
from  EXP1  with  the  zonal  wind  response  from  an  identical 
model  run  using  the  original  Rayleigh  friction  profile  of 
McCormack  [2003],  denoted  as  EXPIRE  in  Table  1.  The 
1 1  -year  variation  in  the  EXP  1  zonal  wind  anomalies  (com¬ 
puted  by  subtracting  the  long-term  monthly  mean)  is  much 
larger  than  in  the  EXPIRE  run,  and  compares  favorably 
with  estimates  of  the  solar  cycle  variability  in  equatorial 
zonal  wind  from  meteorological  analyses  [Crooks  and  Gray, 
2005].  This  result  demonstrates  that  the  increased  Rayleigh 
friction  in  the  earlier  CHEM2D  study  of  McCormack  [2003] 
limited  the  model’s  upper  stratospheric  dynamical  response 
to  solar  UV  variations.  The  implications  of  this  result  for 
modeling  the  solar  cycle  modulation  of  the  QBO  period  are 
discussed  further  in  section  4.2. 

[37]  In  summary,  the  EXP1  modeled  ozone  response  due 
solely  to  changes  in  solar  UV  irradiance  does  not  reproduce 
the  observed  vertical  structure  in  the  1 1 -year  ozone  response 
in  the  stratosphere  at  low  latitudes.  In  the  following  sec¬ 
tions,  we  present  results  from  additional  CHEM2D  experi¬ 
ments  that  include  the  effects  of  the  zonal  wind  QBO  and 
decadal  variability  in  planetary  wave  forcing.  These  results 
show  that  interactions  between  solar  UV  changes,  the  QBO, 
and  planetary  wave  forcing  can  improve  the  agreement 
between  the  modeled  and  observed  11 -year  variation  in 
stratospheric  ozone. 

4.2.  Experiment  2:  Solar  UV  and  QBO  Effects 

[38]  This  section  presents  results  from  a  150-year 
CHEM2D  simulation  that  includes  both  solar  UV  and 
QBO  effects  (EXP2).  This  simulation  demonstrates 
the  solar  cycle  modulation  of  the  zonal  wind  QBO  and 
its  impact  on  stratospheric  ozone  and  temperature.  The 
methodology  of  EXP2  is  similar  to  an  earlier  study  by 
McCormack  [2003]  that  used  a  version  of  CHEM2D  with  a 
2.6  km  vertical  grid  spacing  and  a  model  top  near  ~  106  km. 
This  study  found  a  small  (~  1  month)  decrease  in  the  length 
of  the  QBO  westerly  phase  during  solar  maximum  that 


was  qualitatively  similar  to,  but  much  smaller  than,  the 
3-6  month  decrease  in  the  QBO  westerly  phase  duration  at 
solar  maximum  reported  by  Salby  and  Callaghan  [2000]. 
The  present  version  of  CHEM2D  produces  a  larger  solar 
cycle  modulation  of  the  QBO  westerly  phase  than  that 
reported  by  McCormack  [2003], 

[39]  Figure  7  compares  time  series  of  zonal  winds  over 
the  equator  at  the  32  hPa  level  (~26  km)  from  years  1-51 
of  the  EXP2  (i.e.,  varying  solar  UV,  gray  curve)  and  EXPO 
(i.e.,  fixed  solar  UV,  black  curve)  experiments.  Both 
simulations  are  initialized  with  the  same  model  constituent 
and  transport  fields,  and  they  both  begin  with  the  same 
UV  irradiance  values  corresponding  to  solar  maximum.  We 
find  that  the  two  time  series  start  to  differ  by  year  6, 
corresponding  to  solar  minimum  in  EXP2,  and  remain  out 
of  phase  throughout  this  period  owing  to  differences  in  the 
length  of  the  QBO  westerly  phase  between  solar  maximum 
and  solar  minimum  in  EXP2.  This  behavior  is  in  better 
agreement  with  the  observed  solar  cycle  modulation  of  the 
QBO  than  the  earlier  results  of  McCormack  [2003]  that 
produced  a  weaker  solar  cycle  modulation  of  the  QBO 
westerly  phase. 

[40]  To  quantify  the  solar  cycle  modulation  of  the  model 
zonal  wind  QBO,  we  have  calculated  the  lengths  of  the 
individual  easterly  and  westerly  QBO  phases  at  32  hPa  in 
EXP2  occurring  within  a  20-month  window  centered  on 
solar  maximum  and  minimum.  Figure  8  plots  the  duration 
of  the  individual  easterly  and  westerly  QBO  phases  at 
32  hPa  from  the  entire  150  years  of  EXP2  output  for  solar 
minimum  and  solar  maximum  conditions.  We  find  that  the 
average  westerly  QBO  phase  at  32  hPa  is  15  months  at 
solar  maximum  and  18  months  at  solar  minimum.  A 
Student’s  /-test  shows  that  this  3  month  difference  is 
statistically  significant  at  the  98%  confidence  level.  This 
result  is  in  better  agreement  with  the  observed  3-6  month 
difference  reported  by  Salby  and  Callaghan  [2000]  than  the 
1-month  difference  reported  earlier  by  McCormack  [2003]. 

[41]  The  proposed  origin  of  the  modeled  solar  cycle 
modulation  of  the  westerly  QBO  phase,  which  has  been 
described  previously  by  McCormack  [2003],  is  summarized 
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Figure  8.  Lengths  of  individual  easterly  (E)  and  westerly 
(W)  model  zonal  wind  QBO  phases  at  32  hPa,  sorted 
according  to  the  phase  of  the  solar  cycle,  from  experiment  2. 
The  solar  irradiance  intervals  are  defined  to  be  within 
±20  months  of  the  maximum  or  minimum  in  UV  irradiance. 
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Figure  9.  Annual  average  solar  regression  coefficients  for 
ozone  (percentage  change,  solar  maximum  minus  solar 
minimum)  from  experiment  2,  which  includes  both  solar 
UV  and  QBO  effects.  Shaded  regions  are  not  statistically 
significant  at  the  2<r  level. 

here.  Increased  heating  at  solar  maximum  in  the  upper 
stratosphere  produces  stronger  upward  vertical  motion  near 
the  summer  subtropical  latitudes  and  stronger  meridional 
flow  across  the  equator  from  the  summer  to  winter  hemi¬ 
sphere.  The  anomalous  advection  of  easterly  momentum 
from  the  summer  hemisphere  by  the  meridional  component 
of  the  residual  circulation  offsets  the  westerly  zonal  wind 
tendencies  produced  by  the  parameterized  gravity  wave 
drag  that  drives  the  SAO  westerly  phase.  This  reduces  the 
westerly  zonal  wind  tendency  to  give  overall  weaker 
westerly  flow  in  the  equatorial  upper  stratosphere  during 
solar  maximum  [see  McCormack,  2003,  Figures  6  and  5], 

[42]  Since  the  model  QBO  westerlies  tend  to  emerge 
during  the  SAO  westerly  phase,  solar-induced  changes  in 
the  strength  of  the  upper  stratospheric  westerly  flow  will 
affect  the  timing  and  duration  of  the  QBO  westerly  phase. 
The  model  results  indicate  that  the  weaker  westerly  zonal 
wind  tendencies  under  solar  maximum  conditions  favor  a 
more  rapid  transition  to  the  QBO  easterly  phase,  and  thus  a 
shorter  QBO  west  phase.  In  contrast,  the  timing  and 
duration  of  the  modeled  QBO  easterly  phase  is  affected 
primarily  by  the  annual  cycle  in  tropical  upwelling  in  the 
lower  stratosphere.  This  is  also  seen  observationally  in  the 
preferential  “stalling”  of  the  QBO  easterly  phase  during 
NH  winter.  As  a  result,  the  modeled  dynamical  response  to 
varying  solar  UV  in  the  tropical  upper  stratosphere  affects 
the  duration  of  the  QBO  west  phase  more  than  it  affects  the 
QBO  east  phase. 

[43]  Figures  7  and  8  illustrate  the  solar  cycle  modulation 
of  the  westerly  QBO  phase.  The  high-vertical-resolution 
version  of  CHEM2D  employed  in  the  present  study  produces 
a  larger  average  difference  in  the  length  of  the  QBO 
west  phase  (3  months)  between  solar  maximum  and  mini¬ 
mum  than  was  originally  found  by  McCormack  [2003] 
( 1  month).  The  larger  effect  is  a  consequence  of  the  reduced 
Rayleigh  friction  term  in  the  present  version  of  CHEM2D, 
which  produces  a  larger  solar  cycle  variation  in  model 


dynamics  near  the  stratopause  (Figure  6b)  as  compared  to 
the  earlier  low-vertical-resolution  version  of  the  model. 

[44]  McCormack  [2003]  found  that  the  solar  cycle  mod¬ 
ulation  of  the  QBO  produces  a  quasi-decadal  variation  in 
lower  stratospheric  tropical  upwelling.  This  variation  can 
affect  both  ozone  and  temperature  in  the  tropical  lower 
stratosphere,  where  transport  effects  generally  dominate 
over  photochemical  and  radiative  effects.  The  effect  of  the 
solar-QBO  interaction  on  ozone  can  be  seen  in  Figure  9, 
which  plots  the  annual  average  ozone  solar  regression 
coefficients  derived  from  the  first  50  years  of  EXP2.  The 
presence  of  the  QBO  increases  the  model’s  interannual 
variability,  and  so  the  resulting  2(j  uncertainty  estimates  in 
the  regression  coefficients  are  larger  than  in  EXP1.  Typical 
values  of  the  2er  uncertainty  estimates  for  the  EXP2  ozone 
regression  coefficients  range  from  0.2%  in  the  upper  strato¬ 
sphere  to  0.5%  in  the  lower  stratosphere.  As  with  the  EXP1 
results,  we  find  a  statistically  significant  ozone  response 
throughout  much  of  the  stratosphere  and  lower  mesosphere 
in  the  EXP2  results. 

[45]  Comparing  the  ozone  regression  coefficients  derived 
from  EXP1  (Figure  4)  and  from  the  first  50  years  of  EXP2 
(Figure  9),  we  find  that  the  11 -year  ozone  response  in  the 
tropical  lower  stratosphere  increases  when  the  effects  of  the 
zonal  wind  QBO  are  included.  The  largest  change  is  seen 
between  0  and  10°S  latitude  and  20  and  25  km  altitude, 
where  the  peak  ozone  response  increases  from  (1.5  ± 
0.03)%  in  EXP1  to  (2.2  ±  0.4)%  in  EXP2.  The  fact  that 
this  enhancement  is  not  centered  over  the  equator  is  likely  a 
consequence  of  a  meridional  asymmetry  in  the  modeled 
zonal  wind  QBO  produced  with  the  higher-vertical-resolu¬ 
tion  version  of  CHEM2D.  For  example,  the  peak  zonal 
wind  QBO  amplitude  in  EXP2  is  centered  near  10°S, 
whereas  an  earlier  version  of  CHEM2D  produced  a  zonal 
wind  QBO  centered  over  the  equator  [see,  e.g.,  McCormack 
and  Siskind,  2002,  Figure  4]  that  is  more  consistent  with  the 
observed  latitude  structure  of  the  QBO  [. Baldwin  et  al., 
2001;  Salby  and  Callaghan,  2006]. 

[46]  The  effect  of  the  solar-QBO  interaction  on  temper¬ 
ature  can  be  seen  in  Figure  10,  which  plots  the  annual 
average  temperature  regression  coefficients  derived  from 
the  first  50  years  of  EXP2.  The  EXP2  temperature  response 
is  statistically  significant  at  the  2er  level  throughout  most  of 
the  stratosphere  and  lower  mesosphere,  and  exhibits  a  local 
maximum  in  the  tropical  lower  stratosphere,  similar  to  the 
modeled  ozone  response  in  Figure  9.  The  peak  lower  strato¬ 
spheric  temperature  response  in  EXP2  of  (0.6  ±  0.06)  K  near 
10°S  is  significantly  larger  than  the  corresponding  tempera¬ 
ture  response  of  (0.3  ±  0.01)  in  EXP1  (not  shown). 

[47]  The  presence  of  the  equatorial  zonal  wind  QBO  also 
modifies  the  modeled  solar  cycle  variation  in  extratropical 
dynamics.  For  example.  Figure  1 1  compares  the  composite 
zonal  wind  differences  (i.e.,  solar  maximum  minus  solar 
minimum)  over  the  Northern  Flemisphere  (NH)  extratropics 
in  winter  computed  from  EXP1  (Figure  11a),  and  EXP2 
(Figures  lib- lid).  The  composites  in  Figures  11a  and  lib 
are  constructed  from  50  years  of  model  output  within 
±1  year  of  solar  maximum  or  minimum,  giving  a  total  of 
12  members  in  each  category.  Regions  of  statistically 
significant  differences  are  determined  using  the  Student’s 
/-test.  In  the  case  of  EXP1,  where  only  solar  UV  effects 
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Figure  10.  Annual  average  solar  regression  coefficients 
for  temperature,  expressed  as  the  change  (Kelvin)  between 
solar  maximum  and  solar  minimum,  from  experiment  2  that 
includes  both  solar  UV  and  QBO  effects.  Shaded  regions  are 
not  statistically  significant  at  the  2er  level. 

are  included.  Figure  11a  shows  negligible  zonal  wind 
differences  on  1  December  (Figure  11a,  top).  By  1  January 
(Figure  11a,  middle)  a  2-4  m  s'  1  increase  in  westerly  flow 
at  solar  maximum  develops  in  the  midlatitude  upper  strato¬ 


sphere  and  lower  mesosphere,  which  propagates  poleward 
and  downward  by  1  February  (Figure  11a,  bottom). 

[48]  Similar  behavior  is  noted  in  the  EXP2  zonal  wind 
fields  (Figure  lib),  although  now  the  presence  of  the  zonal 
wind  QBO  introduces  a  higher  degree  of  interannual  vari¬ 
ability  that  reduces  the  statistical  significance  of  the  differ¬ 
ences.  Figure  lib  shows  that  the  QBO  produces  larger 
extratropical  zonal  wind  anomalies  in  NF1  winter  than  in  the 
case  with  no  QBO  (EXPO,  Figure  11a).  The  general  pattern 
of  positive  upper  stratospheric  zonal  wind  anomalies  at 
solar  maximum  appearing  in  early  Northern  winter  and 
propagating  poleward  and  downward  in  the  following 
months  is  similar  to  the  observed  zonal  wind  response 
[Kuroda  and  Kodera,  2002;  Matthes  et  al.,  2004]. 

[49]  To  further  investigate  the  combined  solar  cycle-QBO 
effect  on  the  zonal  winds,  we  sort  the  composite  zonal  wind 
differences  from  EXP2  according  to  the  sign  of  the  model 
zonal  wind  over  the  equator  at  32  hPa  (~26  km),  using  the 
model  output  from  the  month  closest  to  solar  maximum/ 
minimum  with  the  appropriate  QBO  phase.  This  produces  a 
subset  of  4  years  for  each  case.  We  find  a  more  robust 
extratropical  zonal  wind  response  during  Northern  winter  in 
QBO  east  years  (Figure  11c)  than  in  QBO  west  years 
(Figure  lid). 

[50]  An  examination  of  zonal  wind  anomalies  and 
Eliassen-Palm  flux  calculations  from  EXPO  (i.e.,  QBO  with 
no  solar  effect,  not  shown)  indicates  that  the  wintertime  NH 
stratosphere  exhibits  reduced  equatorward  propagation  of 
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Figure  11.  Composite  zonal  wind  differences  (in  m  s-1,  solar  maximum  minus  solar  minimum) 
computed  on  the  first  day  of  December,  January,  and  February  from  (a)  EXP1  with  solar  cycle  and  no 
QBO,  (b)  EXP2  with  solar  cycle  and  QBO,  (c)  EXP2  for  QBO  east  years  only,  and  (d)  EXP2  for  QBO 
west  years  only.  Contours  are  drawn  every  ±1,  ±2,  ±4,  ±6,  ....  Dashed  contours  indicate  negative  values. 
Light  and  dark  gray  shading  denotes  statistical  significance  at  the  90%  and  95%  confidence  levels, 
respectively. 
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Figure  12.  (a)  Latitude  profile  of  planetary  wave 
1  amplitude  imposed  at  the  model  lower  boundary  for 
conditions  at  solar  maximum  (dashed  line)  and  solar 
minimum  (solid  line)  in  experiment  3.  (b)  Time  series  of 
the  vertical  component  of  the  Eliassen-Palm  flux  at  100  hPa 
for  January,  averaged  between  30  and  90°N  latitude,  from 
experiment  3. 


(a)  CHEM2D  PW  1  Amplitude 


planetary  wave  activity  during  easterly  QBO  phase  [e.g., 
Holton  and  Tan ,  1980].  This  change  in  planetary  wave 
propagation  produces  greater  drag  on  the  zonal  winds 
poleward  of  60°N  and  a  reduction  in  the  planetary  wave 
drag  on  the  zonal  winds  at  midlatitudes.  In  this  way,  the 
model  zonal  wind  anomalies  tend  to  reinforce  the  “pure” 
zonal  wind  response  to  the  11 -year  solar  UV  variations 
(Figure  11a)  during  QBO  east  years,  and  tend  to  oppose  the 
pure  solar  response  during  QBO  west  years.  This  behavior 
is  qualitatively  consistent  with  earlier  observational  studies 
that  found  a  significant  correlation  between  solar  activity 
and  stratospheric  temperatures  or  geopotential  heights  in 
Northern  winter  in  years  when  the  equatorial  zonal  wind 
QBO  was  in  its  easterly  phase  [Labitzke  and  van  Loon, 
1988;  Labitzke,  2005], 

4.3.  Experiments  3  and  4:  Effects  of  Imposed  Planetary 
Wave  Forcing 

[51]  Observational  and  3D  modeling  studies  indicate  that 
the  combined  effects  of  the  solar  cycle  and  the  QBO  on 
stratospheric  winds  in  the  Northern  winter  extratropics  can 
reduce  planetary  wave  activity  at  solar  maximum  and  alter 
the  mean  meridional  circulation  of  the  middle  atmosphere 
[Hood  and  Soukharev,  2003;  Matthes  et  al,  2004;  Gray  et 
al.,  2004;  Palmer  and  Gray,  2005],  Since  planetary  wave 
activity  drives  the  upward  component  of  the  mean  merid¬ 
ional  circulation  in  the  tropical  stratosphere  [Holton  et  al., 
1995],  it  has  been  suggested  that  a  decadal  modulation  in 
planetary  wave  activity  associated  with  solar  activity  could 
contribute  to  the  observed  11 -year  signal  in  lower  strato¬ 
spheric  ozone  seen  in  Figure  5  [Hood,  2004].  The  key  to 
this  proposed  mechanism  is  understanding  the  nonlinear 
interaction  between  the  QBO,  planetary  waves,  and  the 
mean  meridional  circulation  [Salby  and  Callaghan,  2006], 

[52]  Unlike  3D  models  where  the  tropospheric  planetary 
wave  forcing  can  vary  in  response  to  changes  in  the  back¬ 
ground  zonal  winds,  CHEM2D  employs  fixed  planetary 


wave  1  amplitude  at  the  model  lower  boundary  (1.3  km). 
Consequently,  the  results  of  EXP1  and  EXP2  could  be 
missing  a  feedback  effect  where  changes  in  planetary  wave 
forcing  in  response  to  solar-induced  zonal  wind  anomalies 
(e.g.,  Figure  11)  alter  tropical  upwelling  and  ultimately 
produce  an  enhanced  dynamical  response  in  ozone  and 
temperature  in  the  tropical  lower  stratosphere.  To  evaluate 
this  effect,  we  have  performed  a  third  experiment  (EXP3) 
where,  in  addition  to  the  1 1  -year  solar  UV  variation  and  the 
QBO,  we  impose  an  11 -year  variation  in  planetary  wave  1 
amplitude  at  the  model  lower  boundary. 

[53]  The  amplitude  of  the  imposed  planetary  wave  forc¬ 
ing  is  based  on  the  study  of  Hood  and  Soukharev  [2003], 
who  found  evidence  of  a  quasi-decadal  variation  in  NH  eddy 
heat  flux  of  ~  1 0%  that  was  anticorrelated  with  the  1 1  -year 
solar  cycle.  Hood  and  Soukharev  [2003]  also  found 
indications  of  a  similar,  albeit  smaller,  quasi-decadal 
variation  in  the  Southern  Flemisphere  (SH)  eddy  heat  flux; 
this  result  is  considered  to  be  less  reliable  than  the  NH  eddy 
heat  flux  variation  because  of  the  poor  data  coverage  at  high 
southern  latitudes. 

[54]  These  variations  in  eddy  heat  flux  may  be  an  indirect 
response  of  the  tropospheric  circulation  to  the  11 -year  solar 
cycle,  as  has  been  suggested,  or  they  may  be  the  result  of 
internal  dynamical  variability  unrelated  to  the  solar  cycle.  In 
either  case,  the  goal  of  this  experiment  is  to  quantify  the 
combined  effect  of  the  solar  UV,  QBO,  and  planetary  wave 
variations  on  stratospheric  ozone  and  temperature  within  the 
2D  modeling  framework.  The  results  of  this  experiment 
provide  further  insight  into  the  possible  origin  of  the 
observed  altitude  dependence  of  the  1 1 -year  ozone  response 
and  can  provide  guidance  for  more  detailed  3D  model 
experiments  in  the  future. 

[55]  Figure  12a  plots  planetary  wave  1  amplitudes  used  in 
the  CHEM2D  model  as  a  function  of  latitude.  In  EXP3,  we 
impose  a  10%  variation  in  the  wave  1  amplitude  over  the 
winter  season  in  both  hemispheres  such  that  the  wave  forcing 
is  higher  during  solar  minimum  and  lower  during  solar 
maximum.  (Note:  CHEM2D  experiments  with  fixed  plane¬ 
tary  wave  1  amplitude  use  the  solar  minimum  values  in 
Figure  12a).  Although  the  SH  eddy  heat  flux  variations 
reported  by  Hood  and  Soukharev  [2003]  may  be  less  reliable 
than  the  NH  estimates,  we  apply  the  wave  1  amplitude  to  both 
hemispheres  to  establish  a  theoretical  upper  limit  on  the 
resulting  ozone  and  temperature  variability. 

[56]  To  demonstrate  that  the  imposed  10%  variation  in 
planetary  wave  1  amplitude  is  consistent  with  the  decadal 
variations  in  NH  eddy  heat  flux  reported  by  Hood  and 
Soukharev  [2003],  Figure  12b  plots  a  time  series  of  the  NH 
average  vertical  component  of  the  EP  flux  (Fz)  at  100  hPa 
for  January  from  experiment  3.  The  eddy  heat  flux  is 
proportional  to  Fz  [Andrews  et  al.,  1987],  and  both  are 
commonly  used  to  diagnose  the  strength  of  the  planetary 
wave  forcing  on  the  stratospheric  mean  meridional  circula¬ 
tion.  The  imposed  11 -year  variation  in  planetary  wave 
1  amplitude  produces  a  ~10%  variation  in  Fz  in  NH  winter, 
with  stronger  forcing  at  solar  minimum  and  weaker  forcing 
at  solar  maximum.  This  time  dependence  agrees  qualita¬ 
tively  with  the  3D  modeling  study  of  Matthes  et  al.  [2004] 
that  reported  a  weaker  meridional  circulation  at  solar 
maximum  in  Northern  winter. 
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Figure  13.  Annual  average  solar  regression  coefficients 
for  ozone  (percentage  change  between  solar  maximum  and 
solar  minimum)  from  experiment  3,  which  includes  varying 
solar  UV,  the  QBO,  and  an  imposed  11 -year  variation  in 
planetary  wave  1  amplitude.  Shaded  regions  are  not 
statistically  significant  at  the  2er  level. 

[57]  Figure  13  plots  the  annual  average  solar  regression 
coefficients  for  ozone  from  EXP3.  Overall,  the  latitude  and 
altitude  dependences  of  the  EXP3  regression  coefficients 
are  similar  to  the  EXP2  results,  with  a  peak  ozone  response 
between  35  and  40  km  and  a  secondary  maximum  in  the 
lower  stratosphere  centered  at  10°S.  However,  the  addition 
of  the  planetary  wave  1  forcing  in  EXP3  produces  a  larger 
peak  ozone  response  in  the  tropical  lower  stratosphere  of 
(2.6  ±  0.5)%  than  seen  in  the  EXP2  results  (Figure  9).  The 
larger  EXP3  ozone  response  in  the  tropical  lower  strato¬ 
sphere  is  accompanied  by  a  small  region  of  reduced 
response  (1.5  ±  0.4%)  at  30  km.  We  also  note  a  slightly 
larger  ozone  response  in  the  equatorial  upper  stratosphere. 
Although  the  differences  between  the  lower  stratospheric 
ozone  responses  in  EXP2  and  EXP3  are  not  statistically 
significant,  the  imposed  planetary  wave  forcing  in  EXP3 
produces  better  qualitative  agreement  with  the  observed 
lower  stratospheric  ozone  response.  The  imposed  planetary 
wave  forcing  also  produces  statistically  significant  ozone 
responses  of  2-3%  in  the  NH  between  40  and  70°N 
peaking  near  18  km  and  in  the  Southern  Hemisphere 
poleward  of  80°S  between  20  and  30  km. 

[58]  Figure  14  plots  the  solar  regression  coefficients 
derived  from  the  CHEM2D  model  temperatures  for  EXP3. 
As  in  the  case  of  EXP2,  we  find  that  the  enhanced  lower 
stratospheric  ozone  response  in  EXP3  (Figure  13)  near  10°S 
is  accompanied  by  an  enhanced  lower  stratospheric  temper¬ 
ature  response  of  (0.7  ±  0.4)  K.  Near  the  equatorial 
stratopause,  we  find  that  the  EXP3  temperature  response 
is  somewhat  diminished  compared  to  the  EXP2  result 
(Figure  10).  Most  differences  in  the  temperature  response 
between  EXP2  and  EXP3  are  not  statistically  significant 
from  each  other  below  50  km.  However,  the  EXP3  results  in 
Figure  14  do  show  large  positive  temperature  responses  of 
~  1 .6  K  at  high  latitudes  in  both  hemispheres  between  50  and 
60  km  that  are  significantly  larger  than  the  corresponding 
EXP2  temperature  response. 


[59]  The  larger  modeled  ozone  and  temperature  responses 
in  the  lower  stratosphere  in  EXP3  suggest  that  the  combined 
effects  of  planetary  wave  forcing  and  solar/QBO  effect 
in  EXP3  could  help  explain  the  observed  2-3%  quasi- 
decadal  variability  in  lower  stratospheric  ozone  reported  by 
Soukharev  and  Hood  [2006],  To  determine  the  model’s 
response  to  the  imposed  planetary  wave  variability  in  the 
absence  of  the  QBO,  a  fourth  50-year  simulation  (EXP4) 
was  carried  out  where  only  the  solar  UV  and  planetary  wave 
amplitude  variations  were  included. 

[60]  Figure  15  plots  the  solar  ozone  regression  coeffi¬ 
cients  for  EXP4.  We  find  that  a  peak  lower  stratospheric 
ozone  response  in  EXP4  of  (2.2  ±  0.4)%  that  is  generally 
less  than  the  peak  response  in  EXP3,  which  includes  the 
effects  of  the  QBO.  Although  these  differences  are  not 
statistically  significant  at  the  2 a  level  because  of  the 
relatively  large  amount  of  interannual  variability  in  the 
model  ozone  fields,  we  find  that  including  the  QBO  in 
EXP3  produces  better  agreement  with  the  observed  lower 
stratospheric  ozone  response. 

[61]  Interestingly,  the  equatorial  upper  stratospheric 
ozone  response  in  EXP4  is  slightly  larger  than  in  EXP3, 
in  better  agreement  with  the  observed  upper  stratospheric 
ozone  response.  Regression  analysis  of  the  EXP4  temper¬ 
atures  (Figure  16)  finds  that,  in  the  absence  of  the  QBO,  the 
11 -year  planetary  wave  forcing  produces  a  smaller  upper 
stratospheric  temperature  response  at  low  latitudes  and  a 
larger  upper  stratospheric  temperature  response  at  high 
latitudes  as  compared  to  EXP3.  Regression  analysis  of 
both  the  EXP3  and  EXP4  temperatures  for  December- 
January-February  and  June-July-August  periods  (not  shown) 
indicates  that  these  high-latitude  upper  stratospheric 
temperature  responses  are  largest  in  the  winter  hemispheres, 
similar  to  the  seasonal  dependence  found  in  earlier  analyses 
[. McCormack  and  Hood,  1996;  Hood,  2004], 

[62]  Both  the  EXP3  and  EXP4  temperature  responses 
suggest  that  the  imposed  variation  in  planetary  wave 


Figure  14.  Annual  average  solar  regression  coefficients 
for  temperature  (change  in  Kelvin  between  solar  maximum 
and  solar  minimum)  from  experiment  3,  which  includes 
varying  solar  UV,  the  QBO,  and  an  imposed  11 -year 
variation  in  planetary  wave  1  amplitude.  Shaded  regions  are 
not  statistically  significant  at  the  2 a  level. 
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Figure  15.  Annual  average  solar  regression  coefficients 
for  ozone  (percentage  change  between  solar  maximum  and 
solar  minimum)  from  experiment  4,  which  includes  varying 
solar  UV  and  an  imposed  11 -year  variation  in  planetary 
wave  1  amplitude.  Shaded  regions  are  not  statistically 
significant  at  the  2<r  level. 

1  amplitude  modulates  the  residual  circulation  in  the  upper 
stratosphere,  producing  more  upwelling  at  low  latitudes  and 
more  downwelling  at  high  latitudes  during  solar  maximum. 
This  effect  is  more  pronounced  in  EXP4,  since  the  presence 
of  the  QBO  in  EXP3  introduces  residual  circulation  anoma¬ 
lies  of  alternating  sign  that  act  to  reduce  the  overall  effect  of 
the  planetary  wave  variations  in  the  annual  mean  response. 
The  diminished  equatorial  upper  stratospheric  temperature 
response  in  EXP4  has  consequences  for  ozone  in  the  upper 
stratosphere,  as  we  discuss  further  in  section  5.  In  the 
tropical  lower  stratosphere,  however,  EXP4  produces  a 
smaller  temperature  response  than  in  EXP3,  indicating  the 
importance  of  the  QBO’s  contribution  to  the  quasi-decadal 
dynamical  variability  in  this  region. 

[63]  To  illustrate  this  point,  Figure  17  plots  composite 
annual  mean  differences  (solar  maximum  minus  solar 
minimum)  in  model  vertical  velocity  (  w*)  from  EXP1, 
EXP2,  EXP3,  and  EXP4.  Overall,  the  effects  of  the  plan¬ 
etary  wave  forcing  in  EXP4  produce  the  strongest  upwelling 
in  the  tropical  upper  stratosphere  (above  45  km)  at  solar 
maximum,  and  the  strongest  downwelling  in  the  extratrop- 
ical  upper  stratosphere.  Only  simulations  that  include  the 
QBO,  however,  exhibit  much  of  a  solar  cycle  variation  in 
the  lower  stratosphere  (below  35  km).  Although  statistically 
significant  differences  cannot  be  obtained  because  of  the 
larger  interannual  variability  associated  with  the  QBO,  the 
w*  anomalies  in  EXP2  and  EXP3  suggest  a  decrease  in 
tropical  upwelling  at  solar  maximum  (i.e.,  negative  anomaly 
in  w*)  between  0°  and  20°S  near  25  km.  An  examination  of 
the  seasonal  variation  in  the  w*  anomalies  (not  shown)  finds 
the  strongest  downwelling  at  high  latitudes  in  winter,  which 
is  consistent  with  the  seasonal  dependence  of  the  high- 
latitude  temperature  response  found  in  the  upper  strato¬ 
sphere  in  both  EXP3  and  EXP4. 

[64]  The  following  section  examines  the  altitude  depen¬ 
dence  of  the  modeled  11 -year  ozone  response  in  more 


detail.  We  find  that  changes  in  the  residual  vertical  velocity 
produced  by  the  combined  effects  of  solar  UV  variation,  the 
QBO,  and  quasi-decadal  variability  in  planetary  wave 
forcing  can  in  turn  produce  vertical  structure  in  the  1 1  -year 
stratospheric  ozone  response  similar  to  that  seen  in  the 
SBUV(/2)  record  over  the  1979-2003  period. 

5.  Comparison  With  SBUV(/2)  Ozone  Response 

[65]  As  the  model  results  in  section  4  show,  incorporating 
the  effects  of  the  zonal  wind  QBO  (EXP2)  and  an  additional 
1 1 -year  variation  in  planetary  wave  forcing  (EXP3)  produce 
lower  stratospheric  ozone  responses  between  0°  and  20°S 
that  are  larger  than  the  response  to  solar  UV  variations  alone 
(EXP1).  The  fact  that  the  EXP3  lower  stratospheric  ozone 
response  is  larger  than  the  EXP4  response  indicates  the 
importance  of  the  QBO  in  producing  a  quasi-decadal  signal 
in  lower  stratospheric  ozone. 

[66]  To  provide  a  more  quantitative  assessment  of  the 
ozone  response  among  the  different  experiments,  Figure  18 
compares  the  vertical  profile  of  annual  mean  solar  regres¬ 
sion  coefficients  for  ozone  averaged  between  0°  and  20°S 
from  EXP1  (solid  black  curve),  EXP2  (solid  gray  curve), 
EXP3  (dashed  gray  curve),  and  EXP  4  (dashed  black  curve) 
with  the  corresponding  profile  of  SBUV(/2)  regression 
coefficients  and  their  2cr  uncertainty  estimates  for  the 
1979-2003  period  [Soukharev  and  Hood ,  2006].  Regres¬ 
sion  coefficients  from  EXPO  (black  dotted  curve),  with  no 
solar  UV  variation,  are  also  plotted  in  Figure  18  to  quantify 
the  quasi-decadal  signal  in  low-latitude  ozone  produced  by 
the  interaction  of  the  QBO  with  the  annual  cycle. 

[67]  Comparison  of  the  observed  and  modeled  solar  cycle 
ozone  responses  in  Figure  18  shows  that  the  combined 
effects  of  solar  UV,  the  QBO,  and  planetary  wave  forcing  in 
EXP3  produce  a  peak  lower  stratospheric  ozone  response  of 
~2.5%  between  22  and  24  km  that  is  in  good  agreement 
with  the  observed  lower  stratospheric  response  at  these 


Latitude 

Figure  16.  Annual  average  solar  regression  coefficients 
for  temperature  (change  in  Kelvin  between  solar  maximum 
and  solar  minimum)  from  experiment  4,  which  includes 
varying  solar  UV  and  an  imposed  11 -year  variation  in 
planetary  wave  1  amplitude.  Shaded  regions  are  not 
statistically  significant  at  the  2er  level. 
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Figure  17.  Annual  mean  composite  differences  in  model  vertical  velocity  w*  (mm  s_1  x  10),  solar 
maximum  minus  solar  minimum,  from  (a)  EXP1,  (b)  EXP2,  (c)  EXP3,  and  (d)  EXP4.  The  zero  contour  is 
solid,  positive  contours  are  dotted,  and  negative  contours  are  dashed.  The  ±0.1  contours  are  omitted 
above  45  km.  Light  and  dark  gray  shading  indicates  regions  where  differences  are  statistically  significant 
at  the  90%  and  95%  confidence  levels,  respectively. 


latitudes.  The  EXP4  results  show  the  largest  ozone  response 
in  the  upper  stratosphere  between  45  and  50  km.  None  of 
the  model  experiments  with  imposed  solar  UV  variations 
are  able  to  reproduce  the  negligible  response  seen  in 
SBUV(/2)  ozone  between  30  and  40  km,  although  it  appears 
that  dynamical  variability  related  to  the  QBO  and  planetary 
wave  forcing  (i.e.,  results  from  EXP2  and  EXP3  in  Figure  1 8) 
acts  to  reduce  the  11 -year  ozone  response  to  solar  UV 
changes  alone  (EXPO)  in  this  altitude  region. 

[68]  Figure  1 7  shows  that  the  enhanced  lower  stratospheric 
ozone  (and  temperature)  response  in  EXP3  can  be  associ¬ 
ated  with  a  decrease  in  the  tropical  upwelling  at  solar 
maximum  (i.e.,  negative  anomaly  in  vv*)  between  0°  and 
20°S  near  25  km.  The  reduced  upwelling  means  there  is 
less  ozone-poor  air  transported  up  from  lower  altitudes  and 
less  adiabatic  cooling  in  this  region. 

[69]  Figure  17  also  shows  that  the  larger  upper  strato¬ 
spheric  ozone  response  in  both  EXP3  and  EXP4  is  related  to 
a  positive  anomaly  in  w*  at  solar  maximum  over  the  equator 
near  50  km.  This  anomaly  indicates  increased  upwelling 
and  adiabatic  cooling  that  acts  to  reduce  the  overall  positive 


temperature  response  caused  by  increased  ozone  heating  at 
solar  maximum.  As  ozone  is  photochemically  controlled  in 
this  region,  the  dynamically  induced  reduction  in  tempera¬ 
ture  at  solar  maximum  leads  to  more  ozone  because  of  the 
inverse  relation  between  ozone  and  temperature  in  the  upper 
stratosphere  [Jucks  and  Salawitch,  2000].  This  anomaly 
also  explains  the  reduced  upper  stratospheric  temperature 
response  at  low  latitudes  in  EXP3  and  EXP4  (Figures  14 
and  16,  respectively)  as  compared  to  EXP2  (Figure  10). 

[70]  While  the  improved  agreement  between  the  modeled 
and  observed  11 -year  ozone  response  shown  in  Figure  18  is 
encouraging,  the  overall  vertical  and  meridional  structure  of 
the  modeled  ozone  response  still  differs  significantly  from 
the  SBUV(/2)  response  in  other  areas,  most  notably  with 
regard  to  the  near  zero  ozone  response  that  is  observed  near 
35  km.  Lee  and  Smith  [2003]  noted  that  this  feature  could 
be  an  artifact  of  volcanic  effects  in  the  observational  record. 
Although  our  experiments  do  not  include  the  impact  of 
volcanic  aerosols  on  ozone,  we  can  conclude  that  dynamical 
variability  related  to  the  QBO  and  imposed  changes  in 
planetary  wave  1  amplitude  is  not  likely  to  reduce  the 


14  of  17 


D16109 


MCCORMACK  ET  AL.:  SOLAR-QBO  IMPACT  ON  STRATOSPHERIC  OZONE 


D16109 


A03  (Solar  Max  -  Solar  Min) 


Figure  18.  Annual  mean  solar  regression  coefficients 
from  model  experiments  0-4  averaged  between  0°  and 
20°S,  and  corresponding  solar  regression  coefficients 
(diamonds)  derived  from  the  1979-2003  SBUV(/2)  record 
with  ±2cr  uncertainty  estimates  [Soukharev  and  Hood, 
2006], 

“pure”  ozone  response  to  solar  UV  changes  near  35  km  by 
more  than  0.5%. 

6.  Summary  and  Discussion 

[71]  Results  from  a  150-year  CHEM2D  model  simulation 
show  that  the  average  duration  of  the  westerly  QBO  phase 
is  3  months  shorter  at  solar  maximum  than  at  solar  mini¬ 
mum.  We  find  a  similar  effect  in  a  50-year  model  simulation 
that  includes  an  imposed  variation  in  planetary  wave  forcing 
in  addition  to  the  solar  UV  and  QBO  effects.  These  results 
confirm  the  initial  findings  of  McCormack  [2003]  that  the 
1 1  -year  cycle  in  solar  UV  irradiance  can  influence  the  period 
of  QBO  in  equatorial  zonal  winds.  We  find  that  the  solar- 
QBO  interaction  is  more  robust  in  the  present  version  of 
CHEM2D  with  higher  vertical  resolution  compared  to  the 
McCormack  [2003]  results  because  of  the  reduced  influence 
of  the  Rayleigh  friction  term.  However,  the  modeled  solar 
cycle  modulation  of  the  westerly  QBO  period  is  still  smaller 
than  the  observed  effect.  Future  studies  should  investigate 
whether  or  not  additional  feedback  processes  involving 
changes  in  parameterized  equatorial  wave  dissipation  due 
to  ozone  heating  anomalies  [Cordero  and  Nathan,  2005]  are 
capable  of  amplifying  the  modeled  solar-QBO  interaction. 

[72]  We  also  quantify  the  response  of  stratospheric  ozone 
and  temperature  to  changes  in  solar  UV  irradiance  by 
applying  a  multiple  linear  regression  model  to  50  years  of 
CHEM2D  output.  First,  considering  only  solar  UV  varia¬ 
tions,  we  reproduce  the  previously  reported  discrepancies 
between  the  modeled  and  observed  altitude  dependence  of 
the  stratospheric  ozone  response.  Specifically,  the  model 
ozone  response  peaks  between  35  and  40  km,  while  the 
observed  response  exhibits  two  peaks  in  the  upper  strato¬ 
sphere  (40-50  km)  and  lower  stratosphere  (below  25  km) 
[e.g.,  Shindell  et  al.,  1999;  Hood,  2004].  Adding  the  effects 
of  an  interactive  parameterization  for  the  zonal  wind  QBO 


produces  a  modeled  lower  stratospheric  ozone  response  of 
~2%  that  improves  the  agreement  with  the  observed  low 
latitude  ozone  response  [ Soukharev  and  Hood,  2006]. 

[73]  Including  the  effects  of  the  QBO  in  the  model 
produces  a  similar  enhancement  in  the  lower  stratospheric 
temperature  response  of  ~0.4  K  as  well  as  enhanced  zonal 
wind  anomalies  (i.e.,  solar  maximum  minus  solar  minimum) 
in  the  Northern  winter  extratropical  stratosphere.  Sorting 
these  zonal  wind  anomalies  according  to  the  phase  of  the 
QBO  shows  that  the  poleward  and  downward  propagation 
of  westerly  anomalies  in  the  tropical  upper  stratosphere 
throughout  the  Northern  winter  during  solar  maximum  is 
most  pronounced  during  the  QBO  east  phase.  This  result 
supports  earlier  observational  studies  that  found  a  higher 
correlation  between  solar  activity  and  polar  stratospheric 
temperatures  during  QBO  east  years  [e.g.,  Labitzke,  2005], 

[74]  We  also  consider  possible  feedback  between  solar- 
induced  changes  in  background  zonal  mean  flow  and 
planetary  wave  propagation  in  the  winter  extratropics  by 
combining  the  effects  of  solar  UV  variations  and  the  zonal 
wind  QBO  with  an  imposed  1 1  -year  variation  in  planetary 
wave  1  amplitude  based  on  the  observational  estimates  of 
Hood  and  Soukharev  [2003],  The  combination  of  smaller 
planetary  wave  forcing  at  solar  maximum  with  the  solar  UV 
variations  and  the  QBO  produces  further  enhancement  of 
the  lower  stratospheric  ozone  response.  It  also  slightly 
reduces  the  ozone  response  between  30  and  40  km  and 
increases  the  ozone  response  between  45  and  55  km.  This 
experiment  produces  the  best  overall  agreement  with  the 
observed  altitude  dependence  of  the  low-latitude  ozone 
response  (Figure  5).  The  imposed  planetary  wave  forcing 
does  not  improve  the  modeled  ozone  response  in  the 
absence  of  the  QBO.  One  remaining  discrepancy  is  that 
none  of  the  experiments  are  able  to  reproduce  the  near  zero 
ozone  response  found  in  the  SBUV(/2)  record  between  30 
and  40  km. 

[75]  The  improved  agreement  between  the  altitude 
dependence  of  the  modeled  and  observed  11 -year  ozone 
response  at  low  latitudes  results  from  changes  in  strato¬ 
spheric  upwelling  throughout  the  stratosphere  (Figure  17). 
In  the  lower  stratosphere,  where  ozone  is  controlled  mainly 
by  transport,  reduced  upwelling  at  low  latitudes  during  solar 
maximum  produces  an  enhancement  in  both  ozone  and 
temperature.  In  the  upper  stratosphere,  where  ozone  is 
photochemically  controlled,  increased  upwelling  at  low 
latitudes  during  solar  maximum  leads  to  a  reduced  positive 
temperature  response  (through  adiabatic  cooling)  and  thus 
an  enhanced  positive  ozone  response. 

[76]  Although  the  2D  model  experiments  presented  here 
do  not  exactly  reproduce  the  altitude  and  latitude  structure 
in  the  observed  11 -year  stratospheric  ozone  response,  our 
findings  indicate  that  a  model  with  interactions  between 
solar  UV  variations,  the  zonal  wind  QBO,  and  changes  in 
planetary  wave  forcing  of  the  middle  atmosphere  circulation 
can  simulate  key  features  of  the  SBUV(/2)  ozone  response 
that  have  been  previously  unexplained.  It  should  be  empha¬ 
sized  that,  unlike  3D  general  circulation  models,  CHEM2D 
cannot  simulate  wave-mean  flow  interactions  in  a  fully  self- 
consistent  manner.  While  the  imposed  11 -year  variation  in 
planetary  wave  1  amplitude  used  in  EXP3  and  EXP4  is 
consistent  with  the  hypothesis  that  changes  in  solar  UV  can 
produce  a  weaker  planetary  wave  forcing  of  the  residual 


15  of  17 


D16109 


MCCORMACK  ET  AL.:  SOLAR-QBO  IMPACT  ON  STRATOSPHERIC  OZONE 


D16109 


circulation  at  solar  maximum,  a  nonsolar  origin  for  the 
planetary  wave  forcing  is  also  possible. 

[77]  Ultimately,  a  3D  modeling  approach  that  includes  an 
internally  generated  QBO  and  SAO  in  equatorial  strato¬ 
spheric  zonal  winds  [e.g.,  Giorgetta  et  al.,  2006]  would  be 
ideal  for  investigating  possible  dynamical  feedbacks  that 
can  contribute  to  the  1 1  -year  ozone  response.  As  the  results 
of  this  study  indicate,  it  is  likely  that  current  discrepancies 
between  modeled  and  observed  solar  cycle  variations  in 
stratospheric  ozone  can  be  reconciled  as  we  improve  our 
understanding  of  the  decadal  variability  in  stratospheric 
dynamics  produced  by  coupling  between  the  solar  cycle, 
the  QBO,  and  planetary  wave  activity. 
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